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Abstract 

All orientation preserving isometries of the hyperbolic three-space are studied, 
and the probability density of conjugate pair separations for each isometry is pre- 
sented. The study is relevant for the cosmic crystallography, and is the theoretical 
counterpart of the mean histograms arising from computer simulations of the isome- 
tries. 



1 Introduction 

Cosmic crystallography is a method to help finding the geometry and the topology of 
the universe fl]. In a close analysis of the method, a description has been presented of 
how each isometry of the universe gives its individual contribution to a pair separation 
histogram (PSH) 0. More recently, the isometries of the infinite 3D euclidean space were 
studied in some detail, and the expected (theoretical) individual contribution of each 
isometry to the PSH was described 0. 

The present report investigates the orientation preserving isometries of H 3 , the 3D 
infinite hyperbolic space with positive definite metric and unitary radius of negative cur- 
vature. To generalize for arbitrary negative curvature one simply needs dividing every 
quantity with dimension length by the radius of curvature. All the results obtained clearly 
reproduce their euclidean counterparts when the radius of curvature tends to infinity. 



2 The probability density of conjugate pair separa- 
tions 

In H 3 we assume a spherical solid ball B. Under an isometry g of H 3 the ball occupies 
a new position B g \ we only consider isometries such that the balls B and B g intersect. 
Assuming a point P G B and denoting as P g G B g its (^-transported, we call the pair 
(P, P g ) a g-pair. We focus our attention on the g-pairs such that P g G BnB g , and assume 
an infinity of points P g uniformly distributed throughout the intersection B fl B g . For the 
given isometry g, we ask for the probability V g (l)dl that a randomly selected g-pair has 
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hyperbolic separation lying between the values I and I + dl; the probability density Vg(l) 
clearly satisfies the normalization condition 



r2a 



where 2a is the diametre of the balls. 



3 Some basic formulas 

A few useful formulas of the hyperbolic trigonometry in 2D are worth having at hand. In a 
geodetic triangle with sides measuring a, b, c and corresponding opposite angles measuring 
a, j3, 7, we have 



the law of sines 



sin a sin f3 sin 7 



sinh a sinh b sinh c ' 
the first law of cosines (cyclic) 

cosh a = cosh b cosh c — sinh b sinh c cos a; 

the second law of cosines (cyclic) 

cos a = — (cos /3 cos 7 — sin (3 sin 7 cosh a) . 



(3) 



(4) 



We also need a few 2D relations between lengths of arcs of geodesies, horocycles and 
equidistant curves to a geodesic; see figure 1 for visualization: 



e' = gcoshr, h' = /icoshr, h = 2 sinh (g/2), 



cosh(#/2), 



tanhp' = tanhpeoshg, 



and as a consequence 



sinh (<//2) = sinh (g/2) cosh r. 



(5) 



(6) 



(a) 



P' 



g 

Figure 1 Some geometric objects in the plane hyperbolic geometry. 

(a) g,g', and r are geodetic arcs; h and /i' are arcs of horocycles; e' is an arc equidistant to g; 
the arcs r are perpendicular to g and to e'; the geodetic arc / is orthogonal to both g and h. 

(b) a geodetic quadrilateral with three right angles, assuming sinh j> sinh g < 1. 
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We often use the line element of H 3 in the cylindrical coordinates, 

dl 2 = dp 2 + sinh 2 p d(f) 2 + cosh 2 p d( 2 ; (7) 

we also use 

dl 2 = -dW 2 + dX 2 + dY 2 + dZ 2 , W = Vl + X 2 + Y 2 + Z 2 , (8) 

where 

W = cosh p cosh (, X = sinh p cos 0, F = sinh p sin 0, Z = cosh p sinh ( . (9) 

In these coordinates the separation / between two points Pi, P2 is given by 

cosh/ = WiW 2 — {XiX 2 + Y\Y 2 + Z 1 Z 2 ) 

= cosh pi cosh p 2 cosh(^! — — sinh P\ sinh p 2 cos(0i — 2 ) • (10) 

For future reference, consider the special situation of the three points (see figure 2) 

Pi =(p 1; 0,0)= (cosh pi; sinh pi, 0,0) , 
Q =(0,0,9)= (cosily 0,0, sinh #) , 

P2 = (P2, 0, g) = (cosh p 2 cosh g; sinh p 2 cos (p, sinh p 2 sin 0) ; 
the separation / between the points Pi and P 2 is then given by (JlQ|), namely 

cosh / = cosh pi cosh p 2 cosh g — sinh pi sinh p 2 cos . (11) 




Figure 2 Relative position of points Pi and P 2 with hyperbolic separation I. 
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4 The isometries of H 3 

We are presently interested in the isometries of H 3 that preserve the orientation. These 
isometries can be classified as 

• hyperbolic, with 5 parametres; these isometries bear some similarity with the 3- 
parametric euclidean translations; 

• elliptic, also with 5 parametres; they are analogous to the also 5-parametric eu- 
clidean rotations in ordinary space; 

• screw motions, with 6 parametres; bear some similarity with the also 6-parametric 
euclidean screw motions; 

• parabolic, with only 4 parametres; again they remind us of the euclidean transla- 
tions. 



In each of these isometries the intersection of the balls B and B g is a rotationally symmetric 
solid lens, whose thickness T, diametre D = 2R, and volume V® we now seek. We denote 
as a the radius of the balls, and m = 2M the separation between their centres C and C g ; 
then we clearly have M < a and (see figure 3) 

T = 2(a-M). (12) 

The separation m depends on the isometry g one is concerned with. Noting that a, M, 
and R make a right-angled triangle with hypotenuse a, we find 

cosh a = cosh M cosh R . (13) 




Figure 3 The balls B and B g with radius a intersect in a solid lens with equatorial radius R 
and thickness T. The centres of the balls are separated m = 2M. 

To have the volume V g B of the solid lens B H B g we first consider a compact cylindrical 
surface C y embeded in the lens, and whose axis coincides with that of the lens (see figure 
4); all points of C y are at a fixed distance y from the axis, so the geometry on C y is 2D- 
euclidean. Denoting as e the length of the generatrices (arcs equidistant to the geodetic 
axis), the area of C y is 

S(y) =2nesmhy. (14) 
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Figure 4 Sketch of a compact cylinder C y inscribed in the solid lens BDB g ; it has radius y and 
generatrices e. 

We must now relate e with a, M, and the variable radius y. In figure 4 we note a right 
angled triangle with sides a (hypotenuse), y, and M + x, so we have 



cosh a = cosh y cosh(M + x) ; 
since from eq. (5a) we have e/2 = x cosh y , then e = 2 cosh y 



cosh" 1 f ss%2 

\ cosh y 



S(y) = 27rsinh2?/ 
The volume V of the lens is clearly 



, i / cosh a \ 
cosh" 1 — — - M 
\ cosh y J 



M 



(15) 
, and 

(16) 



R 



2tt 



S(y)dy 

(tanh a — tanh M) cosh 2 a — (a — M) 



(17) 



It can be checked that when M = we get the hyperbolic volume 7r(sinh2a — 2a) of 
a solid ball with radius a, as expected. Also note that for small values of a and M we 
recover the euclidean volume (27r/3)(a — M) 2 (2a + M) of the solid lens [0] . 



5 Special translations 



Preceding the study of the general hyperbolic isometry of if 3 we first consider the very 
special situation in which the axis L of the isometry g crosses the centre C of the ball B. 
With a = the radius of B, and t = the value of the translation along the axis, we assume 
t < 2a to have nonvanishing intersection B H B g . In this special isometry we clearly have 
the equality m = t. 

According to eq.(^) and figure 1, a point P at a distance r from the axis C is displaced 
under t to a distance I given by 



sinh (1/2) = sinh (t/2) cosh r : 



(18) 



this is a relation involving the variable I (displacement of P), the variable r (distance from 
P to the axis £), and the parametre t (the unique relevant one in this special isometry). 

We next introduce the probability Qg(r)dr that a randomly chosen point P g which is 
in both B and B g be in a radial position between r and r + dr. The probability density 
Qq(r) clearly is proportional to the area S®(r) of the cylinder C r inscribed in the solid 
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lens BnBg (see eq.fllSD), the coefficient of proportionality being the inverse of the volume 

i I cosh a \ 



V® of the lens: 



Q B 9 {r) 



2tt , 
— - smh 2r 

9 



cosh 



y cosh r 



(19) 



The equality of the probabilities Vf(l)dl and Q®{r)dr then gives, using r(l) obtained from 
(0), 



;Q B M 



dr 

dl^ 91 
1 7r sinh I 

Vf sinh 2 t/2 



cosh" 



/ coshasinht/2\ t 
V sinhZ/2 J ~ 2 



(20) 



In the figure 5 we have four instances of Vg(l). Each plot starts abruptly on I = t and 
vanishes when sinh(//2) = coshatanh(t/2). They greatly differ from that of an euclidean 
translation, where Vg(l) — 6(1 — t), a Dirac 5. 




Figure 5 Probability densities V B (l) of pair separations for hyperbolic motions t of a ball B 
with radius a, when the axis of the isometry crosses the centre of the ball. In (a) we took t = 0.3 
and a = 1.5; in (6), t = 1 and a = 2; in (c), t = 3 and a = 4; and in (d), t = 8.5 and a = 4.5. 
Plots (a) and (d) loosely resemble a Dirac 5 mainly because the corresponding solid lenses BnB g 
have volume V® small in comparison with the volume 7r(sinh 2a — 2a) of B. All integrated areas 
are unitary. 
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6 Special screw motions 

It is very simple to generalize the probability density (|20|) to further have a rotation uo of 
B around the axis C of the translation. See figure 6. 




Figure 6 A screw motion of H 3 with axis along the £ axis. Under the preliminary hyperbolic 
motion t (measured on the £ axis) a point P is displaced to the intermediate position P'; then 
the rotation u around £ brings P' to the final position P g geodetically separated I from P. 



The separation / between a point P and the corresponding P g is now given by (|TT 
where we replace p\ = p2 — > r, g — > £, — > c<j: 



cosh Z = cosh r cosh t — sinh r cos u) . 



(21) 



For fixed t and this gives r(Z), from which we derive dr/dl. Since neither the volume 
V g B nor the areas (r) depend on in this special screw motion of B, the density Qg(r) 
is again given by (0). The density Vf(l) is then 



2tt sinh Z 



Vl B (cosht — cosu;) 



cosh 1 coshi 



' cosh t — cos lu \ t 
cosh Z — cos uj I 2 



(22) 



with V g as given in ( |17|) with M = t/2 . A few sample plots of Vz(l) are given in figure 




0.3 



Figure 7 Probability densities Pg(l) of pair separations for special screw motions g = (t,u>) of 
a ball ,8 with radius a, when the axis of the isometry crosses the centre of the ball. In (a) we 
took t = (no translation), a = 1.8, and uj = tt; in (6), £ = 0.3, a = 1.5, and uo = 7r/8; in (c), 
t = 1.2, a = 2.0, and uj = tt/8; and in (d), i = 3.2, a = 2.0, and uj = tt. Each plot starts at 
I = t, and abruptly except when we have pure rotation (t = 0, case (a)). All plots end when 
coshZ = cosh 2 /? cosh t — sinh 2 p cos uj, with cosh/? = cosh a sech(i/2). All integrated areas are 
unitary. 
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We clearly recover the equation ( fffi ) when u = 0. On the other hand, setting t = in 
([22]) gives the P^(0 for a special elliptic isometry, namely a pure rotation u of the ball £> 
when the axis of the rotation contains the centre of the ball: 

vB = Trsinhl cQsh _ 1 / coshasinh(q;/2) 
9 ^ 8 shrV/2) \yWh 2 (//2) -sinhV/2) 

a special P 8 (/) with t = is in figure 7(a). 




7 Parabolic motions 

To describe a parabolic isometry g of if 3 we need first announce its 2-parametric apex 
A, a point at infinity Next we select an arbitrary point C of H 3 , and draw the unique 
horosphere C with centre A and containing C . Then, starting from C we mark an arc of 
horocycle with length /i, laying on C; the direction of the arc and the value of fi demand 
two new parametres and finally fix the isometry g. The horocyclic separation between C 
and its ^-transported C g being /z, the corresponding geodetic separation m is given by 
eq.(5c): 

fi = 2 sinh(m/2) . (24) 

We now consider a solid ball B with radius a and centre C\ clearly there is no loss of 
generality in this last choice. The two parametres (a, m) suffice to completely describe 

vf } {\). 

Denote as r the geodetic altitude of a point P of H 3 relative to the horosphere C; r 
is counted positive if P is outside C, and negative if inside. Also draw the horosphere 
C r with apex A and intersecting P. Under the isometry g all points of C r are equally 
displaced along horocyclic arcs laying on C r , and measuring 

A = /ie r ; (25) 

equivalently, the geodetic separation / between P and P g is given by (see figure 8) 

sinh(//2) = e r sinh(m/2) . (26) 



X 




Cm a 



Figure 8 A parabolic isometry g of H 3 brings the points C and P to C g and P g: respectively; 
m and I are geodetic arcs, and A are horocyclic arcs; r are parallel geodetic arcs orthogonal 
to both /i and A; all arcs lay in a same H 2 . 

For future use we compute alr/dl from (|26|), with m fixed: 

dr 1 , , I . , 
— = -coth - . 27 
all 2 K 2 J K ' 
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The horosphere C r intersects each solid ball B and B g in flat circular disks V r and V rg , 
both with radius p. To have p as a function of r and a we introduce an auxiliary variable 
s (see figure 9) and solve the system 

cosh a = coshscosh(r — /) , e^ = coshs, p = sinhs, (28) 

which gives 

p = J 2(cosh a — cosh r)e r . (29) 




Figure 9 C r and C are equidistant horospheres both intersecting the solid ball B with centre C 
and radius a. We note a right-angled geodetic triangle with sides a (hypotenuse), s, and r — f, 
where / is the geodetic arc orthogonal to both the geodetic arc s and the horocyclic arc p. 



Our geometric situation is now phrased in the following terms: in a 2D flat plane (the 
horosphere C r ) we have two circular disks T> r and T> rg , both with radius p, whose centres 
are separated by A. We ask for the probability Rg(r)dr that a randomly chosen g-pair 
(P, P g ), such that P g G B R B g , has altitude between r and r + dr. Clearly the probability 
density Rg(r) is proportional to the area Sg(r) of the intersection V r r\V rg , the coefficient 
of proportionality being the inverse of the volume V B of the solid lens B R £>„: 



R B g (r) 



y 9 B 



(30) 



The euclidean area Sf (r) is simple to obtain (see figure 10), it is 



S a 9 (r) 



2 ap "4 



yV - A 2 /4 



with A = 2sinh(//2), p(r) as in fl29D , and cosa(r) = A/(2p) . 



(31) 




= P 



Figure 10 D r and P r9 are coplanar flat circular disks with radius p and with separation A 
between their centres C r and C rg . The area Sg(r) of the intersection V r PiV rg is four times the 
shadowed area. 



Since the probabilities Vg(l)dl and Rg(r)dr are the same, we finally have 

m - >,^». , 32 , 
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coth(Z/2) 
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2 cos -i / sinh(//2) \ _ sinh( ^ /2) ^ p2 _ sinh 2 (z/2) 



\ P 



with 



p(l) = ysinh a — [cosh a — sinh(Z/2)/ sinh(m/2)] . 



(33) 



(34) 



See figure 11, where examples of V g (l) for parabolic isometries are reproduced. In each 
plot we have l max and l min given respectively by sinh(//2) = tanh(m/2)e ±R , with coshi? = 
cosh a sech(m/2). 




(c) 



5 / 

Figure 11 Probability densities V g (l) of pair separations for parabolic isometries g of a solid 
ball B with radius a. Under g, the centre of the ball is displaced a geodetic distance m. In (a) 
we took a = 1 and m = 0.5; in (b), a = 8 and m = 15; and in (c), a = 4 and m = 3. All 
integrated areas are unitary. 



8 General translations 

We now generalize the special translations of section 5. We consider a hyperbolic isometry 
g of H 3 whose axis is (, and value t measured along the axis. The solid ball B with radius 
a now has centre C at a distance b from the axis; in section 5 we assumed the special 
value b = 0. Under the isometry g the centre C g of the new ball B g is separated m from 
C; according to (|6]), we have 

sinh(m/2) = sinh(t/2) cosh b. (35) 

We clearly have nonempty intersection B (1 B g only when m < 2a; values of the 
parametres t, b, and a interesting for our purposes then obey the constraint 

sinh(t/2) cosh b < sinha. (36) 

The thickness T, radius R, and volume V g B of the solid lens B PI B g are still given by (|l~2"D, 
(0), and QT7p, with 2M = m(t, 6) as in 

To obtain the probability density P^(/) we follow the same four steps as described in 
ref.|§. The first step is to investigate the shape of the surface B D C r , where C r is the 
infinitely long cylinder with axis ( and radius r. The surface BC\C r is either a topological 
annulus (if 6 + r < a), or a topological disk (if a, b, and r can form a triangle), or is empty 
(if a < \b — r\). To have the dimensions of B fl C r we consider a generic point B = (r, 0, () 
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of its contour; we note that the distance from B to the centre C = (b, 0, 0) of B is the 
radius a, then (|ITD gives C(a, o, r, 0) according to 



cosh a = cosh b cosh r cosh ( — sinh b sinh r cos <fi; (37) 
the variable half width z(<f)) of the intersection is then (eq.(5a) with e' — > z and a — > £) 

z = £ coshr = cosh -1 (a + (3 cos0) coshr, (38) 

where 

cosh a 

a = — — : — , p = tanhotanhr. (39) 

cosh b cosh r 

For fixed values of a, b, and r, the intersection B HC r lies between the curves z((f>) and 
—z{4>). Figure 12(a) depicts an annulus-like intersection B fl C r , which occurs whenever 
< r < a — b (equivalently a > (3 + 1); note that the equator of the annulus measures 
27rsinhr, due to the azimuthal factor = sinh 2 p in (7). Clearly the extremes — ir and 
7r of (f) are identified. 




Figure 12 Sample intersections with a = 2, b = 1, and r = 0.9. 
(a) The annulus-like intersection B (~)C r . 

(6) The (dashed) annulus-like intersection B n £> s n C r when t coshr = 1. 
(c) The (dashed) disk-like intersection B C\ B g H C r when t coshr = 2.5. 

On the other hand, figure 13(a) shows a disk-like intersection B fl C r , which occurs 
whenever a < j3 + 1, with unequal radii z maa; and ma a;Sinhr, with 




Figure 13 Sample intersections with a = 2, 6 = 1, and r = 1.2. 
(a) The disk-like intersection B n C r . 

(6) The disk-like intersection B (~) B g Pi C r when i coshr = 3. 
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The second step is to investigate the shape of the combined intersection B fl B g D C r . 
To examine the possible occurrence of annulus-like intersections BnB g C\C r we project the 
centres of B and B g on the £ axis, and consider the midpoint of these projections; if this 
midpoint lies inside the solid balls, that is, if cosh a > cosh 6 cosh t/2, then annulus-like 
intersections BnB g nC r may occur. Otherwise all intersections B fl B g fl C r are disk-like. 
Clearly B fl B g fl C r is the intersection of B fl C r with i3 5 fl C r , and £> s fl C r is an exact 
copy of BC\C r , only longitudinally displaced a horocyclic distance t coshr along C r . When 
B C\ C r is annulus-like (0 < r < a — b), then B (1 B g C\ C r is either annulus-like (when 
cosh(£/2) < a — (3, see figure 12(6), or disk-like (when a — [3 < cosh(t/2) < a + (3, see 
figure 12(c), or is empty (if cosh(t/2) > a + (3). In the disk-like intersections B fl B g fl C r 
the disk extends from — (p ma x to (p m ax, with (see figures 12(c) and 13(6) ) 



<Pr. 



cos 



_! /cosh(t/2) 



/3 



(41) 



When S fl C r is disk-like, then B nB g nC r is either disk-like (when cosh(t/2) < a + j3, see 
figure 13(6) ), or is empty (if cosh(t/2) > a + (3). 

The third step is to evaluate the area Sg(r) of the surface B fl B g fl C r . To this end we 
define the auxiliary function 

A{a,(3,ip max ) = cosh" (a + (3 cos <f))d<f), (42) 



(43) 



(44) 



in terms of which the areas such as in figures 12(6), 12(c), and 13(6) are 

Sf(r) = [2A(a, (3, ip max ) - ttp max \ sinh2r. 

The fourth and last step is to compute 

v s m _ drjl) S*[r(l)] 
° [) dl Vf ' 

where r(l) is found from (|TS[) , and V g B from (|T7D and (P5]). In figure 14 we reproduce three 
examples of the density V g (l) for general hyperbolic translations; clearly all plots start 
abruptly at l min > t. 

1 







1 




/ 



Figure 14 Sample V g (l) for general translations of H 3 ; all integrated areas are unitary. 

(a) Here a = 3.5, b = 3.1, and t = 2; all intersections BnC r and ,6 n B g n C r are disk-like; 

l mm = t + 0.002. 

(6) Here a = 3.5, b = 2.5, and t = 2; there are disk-like and annulus-like intersections B n C r , 
but all B H B g nC r are disk-like; Z m j n = t. 

(c) Here a = 3.5, 6 = 2.5, and t = 1; there are both types of intersections of both BnC r and 

i3 Pi i3g fl C r , lmin — ^- 
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9 General screw motions 



We already have all elements needed to obtain the probability density Vf(l) of conjugate 
pair separations for a general screw motion, thus generalizing the results of the preceding 
sections 5, 6, and 8. We now make use of all four independent parametres, namely 

• a = radius of the solid balls B and B g , 

• b = distance from the centres of the balls to the axis ( of the isometry, 

• t = translation of the isometry, measured along the axis, and 

• u — angle of rotation of the isometry, around the axis. 

We shall further write all mathematical expressions in a form appropriate for automatic 
calculation of V g (l) in a computer . Without loss of generality for our purposes we assume 
t > and < u < n. 

The separation 2M between the centres of B and B g is now given by ( |T0D 

cosh 2M = cosh 2 b cosh t — sinh 2 b cos u, (45) 

and the condition M < a is necessary to have nonempty intersection B fl B g . Assuming 
this condition is fulfilled, the solid lens B fl B g has thickness T, radius R, and volume V g B 



given by ([12]), (|13"D, and flT7|), respectively, with M given by fl45|). The centre of the lens 
is at a distance a from the axis (, 

tanh6coso;/2 . , 

tanha = — j-*—, (46) 

cosht/2 ' v ; 

and there always exists one diametre of the lens which is directed perpendicular to the axis 
£. The lens intersects the axis whenever a < R, or equivalently cosh6cosht/2 < cosha. 

We next imagine an infinite family of sufficiently long, coaxial (axis (), cylindrical 
surfaces C r with variable radius r. We are interested in the intersection of each C r with 
the solid lens B fl B g \ clearly only values of r in the range (r min , r max ) give nonempty 
intersections B n B g H C r , where 

r m in = (^ - R) 0(<r - R) , r max = a + R. (47) 

Here is the step function with values and 1. Our strategy to approach B fl B g fl C r is 
first study B fl C r , then B g fl C r , and finally the intersection of these two. 

For a given r G (r min , r max ) we note that B fl C r is annulus-like when r < a — b, and 
disk-like otherwise, so we define 

(pmax = 7r&(a~b-r) + cos" 1 f^"^J Q(r + b - a) , (48) 

where a(a, b, r) and /3(b, r) were given in (|39|). Each intersection B fl C r is nonempty for 
4> G (— (f)max, 4>max), and lies between the two curves 

zi = cosh' 1 (a + P cos 4>) cosh rQ(4> + 4> max )Q(4> max - 0), z 2 = -z x \ (49) 

these two curves are drawn on the geometrically flat cylinder C r . 
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The intersection B g D C r is identical to B fl C n but is displaced £ coshr longitudinally 
on C r , and uj azimuthally. It thus lies between the curves 

z 3 = r + cosh^ 1 [a + P cos(0 — uj)] coshr 

X [6(0 -U + (p m ax) ®{<f>max ~ 4>) + ®((frmax ~ <\> + U ~ 2tt) 6(0 + (f> max )), (50) 



where 



2 4 = ~z 3 + 2r, 



r = t coshr 6(0 + max ) 6(0 max - 0); 



(51) 



(52) 



the term with the 6 function containing 2ir in eq. (|50"D is included to allow automatic 
computing. 

As in ref.0, the area S?(r) of the intersection B fl B g fl C r is 



5 g (r) = sinhr 



6(zi — 24) [min(zi, Z3) — max(z 2 , -24)] d0. 



(53) 



Finally, the probability density Vz(l) is given by ( f44|) with r(Z) coming from (21); we find 
that 



dr sinh / csch2r 
dl cosh t — cos uj 

In figure 15 a few sample plots of for screw motions in H 3 are given. 




(54) 



Figure 15 Sample Vg(l) for general screw motions of H ; all integrated areas are unitary, 
(a) Here a = 4, 6 = 3, t = 1.5, and a> = tt/2. 
(6) Here a = 4, 6 = 2.5, i = 1, and u = tt/2. 
(c) Here a = 4, 6 = 2.5, i = 0.5, and u = 7r/4. 



10 Conclusion 



The three plots in figure 15 are the output of a computer program whose inputs are the 
values of a, b, t, and uj; given these inputs, the program proceeds without any intervention. 
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The three parametres (t, u, b) related to the screw motion can be extracted from a 
4x4 matrix Ai g , which expresses the motion in terms of the minkowskian coordinates 
(W; X,Y, Z) Q. Indeed, it can be shown that the trace T, the sum £ of the principal 
minors of order 2, and the time-time coefficient U of the matrix Ai g are 

T = 2(cosht + costu), E = 2(1 + 2 cosh i cos a; ), U = cosh t cosh 2 b — costo sinh 2 6; (55) 



if we define V = y/8 + T 2 — 4E then we obtain 

1 1 4U — T 

cosht = ~(T + V), cosou = ~(T-V), cosh26= — - — . (56) 

To close this report we show through a concrete example how to use the functions 
Vg(l) to get information about the topology of the universe. We need first briefly recall 
the theory that underlies the subject; for details see 0, 0, ||. 

Expected (or theoretical) normalized histograms <p(i) of pair separations are decom- 
posed as 



b un (i) + 



1 



n — 



where <p un (i) is the expected normalized histogram of the uncorrelated pair separations, 
the i denotes an interval of separations (a bin in the histogram), n is the finite number 
of objects in the solid ball B, u g = N g /n = v g -\ with N g the number of g-pairs with both 
members inside B, and 4> 9 (i) is the expected normalized histogram of the g-pairs. These 
4> 9 (i) are the histogramic counterparts of the functions Vg(l) of || and of this report. 
From (|57|), and accepting a suggestion by Fagundes and Gausmann [0], we write 



71 — 1)0 — =171 



£^[0 w -^i+E^-<n (58) 



where <ft sc (i) is the expected normalized histogram of pair separations in a simply con- 
nected ball with same radius and geometry as £>; it is the histogramic counterpart of 
the function jF H (a, Z) of || and of J-b(o>, s) of 0. In the limit n — ► oo the products 
n[(j) — 4> sc ] =: ip B and n[(p un — <f) sc ] =: ip® n remain finite, and we write 

AO = (0 + rf(0, rf(0 : = £ "*P?(0 - ^(01- (59) 

To go from the various histograms f(i) to the corresponding functions f(l) we have 
simply made the number of bins tend to infinity. The function <p B (l) has been called 
the topological signature of a ball B in a multiply connected space; since in practice the 
function (p% n (t) is usually small valued when compared with both (p B (l) and <p®(l), the 
function <pf(l) is generally a good approximation of the topological signature ^ 3 {l). 

We now turn to a specific example: that of a ball B in the Seifert- Weber dodecahedral 
space. This multiply connected hyperbolic three-space is obtained from a regular solid 
dodecahedron T> by pairwise identifying opposite faces using twists of 3/10 of a revolution 



lOfl . We make the centres of B and T> coincide (a rather uncopernic assumption), and 
choose B tangent to the edges of V. 

Computer simulations of (n — 1) [<f)(i) — (f> sc {i)} for the hyperbolic ball B are given in 
the literature (see figure 7 in || or figure 3 in [0]), and we now construct its approximate 
expected counterpart (0> eq- (|59|) ■ 
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We first select two of the 12 matrices of face-pairing isometries of T>: 



Mi 



( 3.736068 

-0.309017 

0.951057 

V 3.599751 



3.599751 ^ 
-0.951057 
-0.309017 

3.736068 J 



(60) 



( 3.736068 -3.219715 -1.609857 ^ 

-3.219715 2.927051 0.425325 1.618034 




0.425325 -0.309017 0.850651 



(61) 



V -1.609857 1.618034 -0.850651 0.5 

sufficient for our purposes. Applying (|56| ) to any of M.i or M. 2 we find the values 

t x = 1.992, ux = 108°, bx = 0, (62) 

for the translation, the rotation, and the distance from the axis of the isometry to the 
centre of B (the origin of coordinates). From ( f45"D and fl62| ) we obtain half separation 
Mx = 0.996 between the centres of B and B g , a value smaller than the radius a = 1.439 
of B. 

We also need consider 60 other isometries, whose common prototype matrix is 



M 3 = MxM 2 



( 8.163119 -6.204554 -3.062131 -4.214661 \ 

0.994947 -0.5 0.162460 -1.309017 

-3.062131 2.915224 0.5 1.275976 

V 7.434376 -5.545085 -3.178089 -3.927051 J 



(63) 



each such isometry gives, from (EB 



t 3 = 1.746, u 3 « 147°, 6 3 = 0.999. (64) 

These 60 isometries also contribute to ([59]), since from ( fiBf) and ( |64|) we find half separation 
M 3 = 1.395, a value smaller than a. All other isometries seem to give M > a, so they do 
not contribute to (|59|) . 

For the 12 fundamental isometries (^) we find intersections B fl B g with volume 
V® = 1.377 as given by (|1~7D ; all produce the same spectrum Vf{l), which has nonzero 
values only in the interval I G (1.99, 2.79). On the other hand, for the 60 isometries ( p4j ) we 
find V® = 0.011433, and a spectrum V®(1) with nonzero values only when I G (1.75, 1.99). 
Finally, the probability density Vf c {l) for pair separations in a hyperbolic ball is the 
function J r H{a, I) of ||, or the function J-h(cl, s) of 0; for unitary curvature of the space, 
the volume of the ball with radius a = 1.439 is Vf c = 18.8. From we then have (see 
figure 16) 



i2v; £ 

sc 



6oy 3 B 

sc 



(65) 
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Figure 16 Approximate topological signature ^p(i) for an observed universe endowed with the 
Seifert- Weber dodecahedral topology and unitary negative curvature. The centre of observation 
and that of the dodecahedron coincide, and the event horizon is supposed a = 1.44 away. The 
discontinuities observed at I = 1.75 and I = 1.99 derive from isometries, as described in the text. 
Their localization and strength are good indicators of the topology of the universe. 



The approximate signature <p®(l) of figure 16 bears close similarity with the corre- 
sponding histograms figure 7 in || and figure 3 in [IT]. However, some small distortion 
can be seen, probably arising because the uncorrelated contribution <p® n (l), present in 
(Ej^), was not taken into account. As a matter of fact, we have not been able to obtain 
the expected </3f n (/) neither for the present Seifert- Weber space nor for any simpler 3D 
nontrivial manifold, such as the three-torus. Even for the two-torus that function has 
been eluding our efforts; only for the one-torus (a circle) we have already succeeded in 
finding the <p B un (l) Q 
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